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The azimuthal anisotropies of the collective transverse flow of charged hadrons are investigated in 
a wide range of heavy-ion collision energies within the microscopic parton-hadron-string dynamics 
(PHSD) transport approach, which incorporates explicit partonic degrees of freedom in terms of 
strongly interacting quasiparticles (quarks and gluons) in line with an equation of state from lattice 
QCD as well as the dynamical hadronization and hadronic collision dynamics in the final reaction 
phase. The experimentally observed increase of the elliptic flow V2 of charged hadrons with collision 
energy is successfully described in terms of the PHSD approach. The PHSD scaling properties of 
various collective observables are confronted with experimental data as well as with hydrodynamic 
predictions. The analysis of higher-order harmonics V3 and V4 in the azimuthal angular distribution 
shows a similar tendency of growing deviations between partonic and purely hadronic models with 
increasing collision energy. This demonstrates that the excitation functions of azimuthal anisotropies 
reflect the increasing role of quark-gluon degrees of freedom in the early phase of relativistic heavy- 
ion collisions. Furthermore, the specific variation of the ratio V4/ (1^2) 2 with respect to bombarding 
energy, centrality, and transverse momentum is found to provide valuable information on the un- 
derlying dynamics. 



PACS numbers: 25.75.-q, 25.75.Ag 
I. INTRODUCTION 

The discovery of large azimuthal anisotropic flow at 
the Relativistic Heavy Ion Collider (RHIC) provides con- 
clusive evidence for the creation of dense partonic matter 
in ultra relativistic nucleus-nucleus collisions. With suffi- 
ciently strong parton interactions, the medium in the col- 
lision zone can be expected to achieve local equilibrium 
and exhibit approximately hydrodynamic flow [1-3] . The 
momentum anisotropy is generated due to pressure gra- 
dients of the initial "almond-shaped" collision zone pro- 
duced in nonccntral collisions [1, 2]. The azimuthal pres- 
sure gradient extinguishes itself soon after the start of 
the hydrodynamic evolution, so the final flow is insensi- 
tive to later stages of the fireball evolution. The pres- 
sure gradients have to be large enough to translate an 
early asymmetry in density of the initial state to a final- 
state momentum-space anisotropy. In these collisions a 
new state of strongly interacting matter is created, being 
characterized by a very low shear viscosity r\ to entropy 
density s ratio, rj/s, close to a nearly perfect fluid [4-6]. 
Lattice QCD (1QCD) calculations [7-9] indicate that a 
crossover region between hadron and quark-gluon mat- 
ter should have been reached in these experiments. 

An experimental manifestation of this collective flow 
is the anisotropic emission of charged particles in the 
plane transverse to the beam direction. This anisotropy 
is described by the different flow parameters defined as 
the proper Fourier coefficients v n of the particle distribu- 
tions in azimuthal angle ip with respect to the reaction 



plane angle ^rp. At the highest RHIC collision energy 
of y'sNN = 200 GeV, differential elliptic flow measure- 
ments V2(pt) have been reported for a broad range of 
ccntralities or number of participants N part . For N part 
estimates, the geometric fluctuations associated with the 
positions of the nucleons in the collision zone serve as the 
underlying origin of the initial eccentricity fluctuations. 
These data are found to be in accord with model calcula- 
tions that an essentially locally equilibrated quark-gluon 
plasma (QGP) has little or no viscosity [2, 10-12]. Collec- 
tive flow continues to play a central role in characterizing 
the transport properties of the strongly interacting mat- 
ter produced in heavy- ion collisions at RHIC. Particle 
anisotropy measurements are considered as key observ- 
ables for a reliable extraction of transport coefficients. 

A quark-number scaling of the elliptic-flow data is ob- 
served for a broad range of particle species, collision cen- 
tralizes, and transverse kinetic energy, which is inter- 
preted as being due to the development of substantial 
collectivity in the partonic phase [13]. Small violations 
of the scaling of V2(N part ) with the initial eccentricity 
of the collision zone €2 suggest a strongly coupled low- 
viscosity plasma 77/s ~ (1 — 2)/(47r) in energetic Au+Au 
collisions [13-15]. The initial eccentricity of the collision 
zone (and its associated fluctuations) has proven to be an 
essential ingredient for these extractions. Nevertheless, 
the degree to which the QGP is thcrmalizcd is still being 
debated [16]. 

It was shown before that higher-order anisotropy har- 
monics, in particular U4, may provide a more sensitive 
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constraint on the magnitude of rj/s and the freeze-out 
dynamics, and the ratio f4/(w2) 2 might indicate whether 
a full local equilibrium is achieved in the QGP [17]. The 
role of fluctuations and so-called nonflow correlations are 
important for such measurements. It is well established 
that initial eccentricity fluctuations significantly influ- 
ence the magnitudes of V2.4 [18]. However, the precise 
role of nonflow correlations, which lead to a systematic 
error in the determination of i>2,4, is less clear. Recently, 
significant attention has been given to the study of the 
influence of initial geometry fluctuations on higher or- 
der eccentricities e„(n > 3) for a better understanding 
of how such fluctuations manifest themselves in the har- 
monic flow correlations characterized by v n . Even more, 
it was proposed that the analysis of for all values of 
n can be considered as an analogous measurement to the 
power spectrum extracted from the cosmic microwave 
background Radiation providing a possibility to observe 
superhorizon fluctuations [19]. 

More recently, the importance of the triangular flow V3, 
which originates from fluctuations in the initial collision 
geometry, has been pointed out [20-22] . The participant 
triangularity characterizes the triangular anisotropy of 
the initial nuclear overlap geometry. It arises from event- 
by-event fluctuations in the participant-nucleon colli- 
sion space-time points and corresponds to a large third 
Fourier component in the two-particle azimuthal correla- 
tions at large pscudo rapidity separation Arj. This fact 
suggests a significant contribution of the triangular flow 
to the ridge phenomenon and broad away-side structures 
observed in the RHIC data [23]. The ridge might be re- 
lated to flux-tube-like structures in the initial state as 
argued in Ref. [24] or successive coherent gluon radiation 
as suggested in Ref. [25]. 

A large number of anisotropic flow measurements 
have been performed by many experimental groups at 
SIS (Schwerionensynchroton) , AGS (Alternating Gradi- 
ent Synchrotron), SPS (Super Proton Synchrotron), and 
RHIC (Relativistic Heavy Ion Collider) energies over the 
last twenty years. Very recently, the azimuthal asymme- 
try has been measured also at the LHC (Large Hadron 
Collider) [26]. However, the fact that these data have 
not been obtained under the same experimental condi- 
tions as at RHIC experiments, does not directly allow 
for a detailed and meaningful comparison in most cases. 
The experimental differences include different centrality 
selection, different transverse momentum acceptance, dif- 
ferent particle species, different rapidity coverage, and 
different methods for flow analysis as pointed out in 
Ref. [27]. 

The Beam Energy Scan (BES) program proposed at 
RHIC [28] covers the energy interval from ^snn = 
200 GeV, where partonic degrees of freedom play a de- 
cisive role, down to the AGS energy of ^/snn ~ 5 GeV, 
where most experimental data may be described success- 
fully in terms of hadronic degrees of freedom, only. Low- 
ering the RHIC collision energy and studying the energy 
dependence of anisotropic flow allows us to search for the 



possible onset of the transition to a phase with partonic 
degrees of freedom at an early stage of the collision, as 
well as possibly to identify the location of the critical end 
point that terminates the crossover transition at small 
quark-chemical potential to a first-order phase transition 
at higher quark-chemical potential [13, 29]. 

This work aims to study excitation functions for differ- 
ent harmonics of the chargcd-particlc anisotropy in the 
azimuthal angle at midrapidity in a wide transient en- 
ergy range, i.e., from the AGS to the top RHIC energy. 
The first attempts to explain the preliminary STAR data 
with respect to the observed increase of the elliptic flow 
vi with the collision energy have failed, since the tradi- 
tional available models did not allow clarification of the 
role of the partonic phase [30]. In this paper, as an ex- 
tension of our recent study in Ref. [31], we investigate the 
energy behavior of different flow coefficients, their scal- 
ing properties and differential distributions. Our anal- 
ysis of the STAR/PHENIX RHIC data - based on re- 
cent results of the BES program - will be performed 
within the parton-hadron-string dynamics (PHSD) trans- 
port model [32] that includes explicit partonic degrees of 
freedom as well as a dynamical hadronization scheme for 
the transition from partonic to hadronic degrees of free- 
dom and vice versa. 

The paper is organized as follows: In Section II we 
will briefly recall the main ingredients of the PHSD ap- 
proach as well as the performance of PHSD for relativistic 
heavy-ion collisions from the lower SPS to the top RHIC 
energies. Section III is devoted to the actual results from 
PHSD for the excitation function of the elliptic flow V2 
in comparison to the hadron-string dynamics (HSD) ap- 
proach and other related models, as well as to the avail- 
able data from the STAR and PHENIX Collaborations. 
We also provide results from the PHSD and HSD mod- 
els for the excitation functions of ^3 and in view of 
the Beam Energy Scan (BES) program at RHIC in order 
to identify partonic contributions. Scaling properties of 
experimental data, in particular the universal and longi- 
tudinal scaling relations found empirically, are elaborated 
here and compared to a hydrodynamic description. To 
be more specific, we will also present the calculated re- 
sults for the pt dependence of elliptic flow at midrapidity 
for minimum bias collisions of Au+Au for ^snn from 5 
to 200 GeV. Furthermore, the centrality dependence of 
i> 2 , V3, and V4 will be addressed at the top RHIC energy. 
Section IV provides the conclusions of our present study 
and indicates the open problems. 



II. THE PHSD MODEL 

The dynamics of partons, hadrons and strings in rel- 
ativistic nucleus-nucleus collisions is analyzed within 
the novel parton-hadron-string dynamics (PHSD) ap- 
proach [32, 33]. In this transport approach the partonic 
dynamics is based on the dynamical quasiparticlc model 
(DQPM) [34, 35], which describes QCD properties in 
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terms of single-particle Green's functions [in the sense 
of a two-particle irreducible (2PI) approach] . In Ref. [33] 
the actual (essentially three) DQPM parameters for the 
temperature-dependent effective coupling have been fit- 
ted to the recent lattice QCD results of Ref. [9]. The 
latter lead to a critical temperature T c sa 160 McV, 
which corresponds to a critical energy density of e c w 
0.5 GeV/fm 3 . In PHSD the parton spectral functions 
Pj (j — 5) 9 ) are no longer ^-functions in the invariant 
mass squared (as in conventional cascade or transport 
models), but are taken as 



(1) 

separately for quarks/antiquarks and gluons (j — q, q,g). 
With the convention i5 2 (p 2 ) = p 2 + M 2 — 7 2 , the pa- 
rameters M 2 and jj are directly related to the real 
and imaginary parts of the retarded self-energy, e.g., 
2i~fjLo. The spectral function (1) is anti- 
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symmetric in u> and is normalized as 

f°° du , . r doj n . . 

/ 7T w Pj( u ' P) = / 7T 2uj Po( u ' P) = 1 • ( 2 ) 

The actual parameters in Eq. (1), i.e., the gluon mass M g 
and width 7 S - employed as input in the PHSD calcula- 
tions - as well as the quark mass M q and width j q , are 
depicted in Fig. I of Ref. [33] as a function of the scaled 
temperature T/T c . As mentioned above, these values 
for the masses and widths have been fixed by fitting the 
lattice QCD results from Ref. [9] in thermodynamic equi- 
librium. We recall that the DQPM allows extraction of a 
potential energy density V p from the space- like part of the 
energy-momentum tensor, which can be tabulated, e.g., 
as a function of the scalar parton density p s . Derivatives 
of V p with respect to p s then define a scalar mean-field 
potential U s {p s ), which enters the equation of motion for 
the dynamical partonic quasiparticles. As one can see 
from Fig. 1, the scalar potential is rather large and non- 
linear ly increases with p s . This implies that the repulsive 
force due to U s (p s ) will change in a nonmonotonous way 
with the scalar density. The vector mean-field potential 
is not negligible too, especially at high p s , and induces 
a Lorentz force for the partons. Note that the vector 
mean-field potential vanishes with decreasing scalar den- 
sity whereas the scalar mean-field potential approaches a 
constant value for p s — s- 0. 

Furthermore, a two-body interaction strength can be 
extracted from the DQPM as well from the quasipar- 
ticle width in line with Ref. [6] (cf. Refs. [32, 33] for 
details). The transition from partonic to hadronic de- 
grees of freedom (and vice versa) is described by co- 
variant transition rates for the fusion of quark-antiquark 
pairs or three quarks (antiquarks), respectively, obey- 
ing flavor current-conservation, color neutrality, as well 
as energy-momentum conservation. Since the dynami- 
cal quarks and antiquarks become very massive close to 
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FIG. 1: (Color online) The scalar and vector mean-field po- 
tentials in the present PHSD model as a function of the scalar 
density p s of partons. 



the phase transition, the formed resonant "pre-hadronic" 
color-dipole states (qq or qqq) are of high invariant mass 
too, and sequentially decay to the ground-state meson 
and baryon octets, increasing the total entropy. 

On the hadronic side PHSD includes explicitly the 
baryon octet and decouplct, the 0~- and I "-meson 
nonets, as well as selected higher resonances as in 
the hadron-string dynamics (HSD) approach [36, 37]. 
Hadrons of higher masses (>1.5 GeV in the case of 
baryons and >1.3 GeV in the case of mesons) are treated 
as "strings" (color dipoles) that decay to the known (low- 
mass) hadrons according to the JETSET algorithm [38] . 
We discard an explicit recapitulation of the string for- 
mation and decay and refer the reader to the original 
work [38]. Note that PHSD and HSD (without explicit 
partonic degrees of freedom) merge at low energy den- 
sity, in particular below the critical energy density s c ~ 
0.5 GeV/fm 3 . 

The PHSD approach has been applied to nucleus- 
nucleus collisions from y/~SN~N ~ 5 to 200 GeV in 
Refs. [32, 33] in order to explore the space-time regions of 
"partonic matter". It was found that even central colli- 
sions at the top SPS energy of ^snn = 17.3 GeV show a 
large fraction of nonpartonic, i.e., hadronic or string- like 
matter, which can be viewed as a hadronic corona. This 
finding implies that neither hadronic nor only partonic 
"models" can be employed to extract physical conclusions 
in comparing model results with data. In addition, we 
have found that the partonic phase has a low impact on 
rapidity distributions of hadrons but a sizable influence 
on the transverse mass distribution of final kaons due 
to the repulsive partonic mean fields [32]. It has been, 
furthermore, demonstrated in Ref. [33] that at ^/s^n = 
200 GeV the PHSD model gives a reasonable reproduc- 
tion of hadron rapidity distributions and transverse mass 
spectra, and also a fair description of the elliptic flow of 
charged hadrons as a function of the centrality of the 
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reaction and the transverse momentum px- 

Furthermore, an approximate quark-number scaling of 
the elliptic flow V2 of identified hadrons is observed in 
the PHSD results at top RHIC energies too. As indi- 
cated above, PHSD merges to HSD in the lower (tran- 
sient) energy regime. Both approaches are well in line 
with experimental data in the lower SPS energy regime 
as shown in Rcf. [32] . All these previous findings provide 
promising perspectives to use PHSD in the whole range 
from about ^/Sjvjv = 5 to 200 GeV. 

III. RESULTS FOR COLLECTIVE FLOWS 
A. Elliptic flow 




is [GeV] 



The largest component, known as elliptic flow V2, is 
one of the early observations at RHIC [39]. More re- 
cently, it was noticed that fluctuations in the initial geom- 
etry are very important [20] . The elliptic flow coefficient 
is a widely used quantity characterizing the azimuthal 
anisotropy of emitted particles, 

2 2 

v 2 =< cos(2iP - 2V RP ) >=< Vj f^\ > , (3) 

Px + Py 

where ^ rp is the azimuth of the reaction plane, p x and p y 
are the x and y component of the particle momenta, and 
the brackets denote averaging over particles and events. 
This coefficient can be considered as a function of cen- 
trality, pscudorapidity 77, and/or transverse momentum 
Pt- We note that the reaction plane in PHSD is given 
by the x — z plane with the z axis in the beam direction. 
The reaction plane is defined as a plane containing the 
beam axes and the impact parameter vector. 

We recall that at high bombarding energies the lon- 
gitudinal size of the Lorentz contracted nuclei becomes 
negligible compared to its transverse size. The forward 
shadowing effect then goes away and the elliptic flow fully 
develops in-plane, leading to a positive value of the av- 
erage flow i'2 since no shadowing from spectators takes 
place. In Fig. 2 the experimental V2 data compilation 
for the transient energy range is compared to the results 
from HSD calculations and further available model re- 
sults as included in Ref. [30]. The centrality selection is 
the same for the data and the various models. 

In order to interpret the results in Fig. 2 we have to 
recall the various ingredients of the models employed 
for comparison. The UrQMD (ultrarclativistic quan- 
tum molecular dynamics) model is a microscopic trans- 
port theory based on the relativistic Boltzmann equa- 
tion [40]. It allows for the on-shell propagation of all 
hadrons along classical trajectories in combination with 
stochastic binary scattering, color string formation, and 
resonance decay. The model incorporates baryon-baryon, 
meson-baryon, and meson-meson interactions based on 
experimental data (when possible). This Boltzmann- like 
hadronic transport model has been employed for proton- 
nucleus and nucleus-nucleus collisions from AGS to RHIC 



FIG. 2: (Color online) The average elliptic flow 1)2 of 
charged particles at midrapidity for minimum bias collisions 
at ^sKw = 9.2, 19.6, 62.4 and 200 GeV (given by stars) is 
taken from the data compilation of Ref. [30] . The correspond- 
ing results from different models are compared to the data and 
explained in more detail in the text. 



energies [40] . The comparison of the data on V2 to those 
from the UrQMD model will thus essentially provide in- 
formation on the contribution from the hadronic phase. 
As seen in Fig. 2, being in agreement with data at the 
lowest energy t/snJj = 9.2 GeV, the UrQMD model re- 
sults then either remain approximately constant or de- 
crease slightly with increasing tJsnn', UrQMD thus does 
not reproduce the rise of V2 with the collision energy as 
seen experimentally. 

The HSD model [37, 41] is also a hadron-string model 
including formally the same processes as UrQMD. How- 
ever, being based on the off-shell generalized transport 
equation [42] followed from Kadanoff-Baym approach, 
the quasiparticles in the HSD model take into account 
in-medium modifications of their properties in the nu- 
clear environment, which is rather essential for many ob- 
scrvablcs and in particular for dileptons. Detailed com- 
parisons between HSD and UrQMD for central Au + Au 
(Pb + Pb) collisions have been reported in Refs. [43, 44] 
from AGS to top SPS energies with respect to a large ex- 
perimental data set. Indeed, both hadronic approaches 
yield similar results on the level of 20%-30%, which is also 
the maximum deviation from the data sets. Accordingly, 
the HSD model also predicts an approximately energy- 
independent flow V2 in quite close agreement with the 
UrQMD results. We may thus conclude that the rise of 
V2 with bombarding energy is not due to hadronic in- 
teractions and models with partonic degrees of freedom 
have to be addressed. 

The AMPT (a multi phase transport) model [45] uses 
initial conditions of a perturbative QCD (pQCD) in- 
spired model, which produces multiple minijet partons 
according to the number of binary initial nucleon-nucleon 
collisions. These (masslcss) minijet partons undergo 
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scattering (without potentials) before they are allowed 
to fragment into hadrons. The string melting (SM) ver- 
sion of the AMPT model (labeled in Fig. 2 as AMPT- 
SM) is based on the idea that the existence of strings 
(or hadrons) is impossible for energy densities beyond 
a critical value of e ~ 1 GeV/fm 3 . Hence they need 
to melt the strings to (massless) partons. This is done 
by converting the mesons to a quark and antiquark pair, 
baryons to three quarks, etc, fulfilling energy-momentum 
conservation. The subsequent scatterings of the quarks 
are based on a parton cascade with (adjustable) effec- 
tive cross sections that are significantly larger than those 
from pQCD [45]. Once the partonic interactions termi- 
nate, the partons hadronize through the mechanism of 
parton coalescence. 

We find from Fig. 2 that the interactions between the 
minijet partons in the AMPT model indeed increase the 
elliptic flow significantly as compared to the hadronic 
models UrQMD and HSD. An additional inclusion of in- 
teractions between partons in the AMPT-SM model gives 
rise to another 20% of t>2, bringing it into agreement (for 
AMPT-SM) with the data at the maximal collision en- 
ergy. So, both versions of the AMPT model indicate the 
importance of partonic contributions to the observed el- 
liptic flow V2 but do not reproduce its growth with ^/sJvw. 
The authors address this result to the partonic equation 
of state (EoS) employed, which corresponds to a massless 
and noninteracting relativistic gas of particles. This EoS 
deviates severely from the results of lattice QCD calcu- 
lations for temperatures below 2-3 T c . Accordingly, the 
degrees of freedom are propagated without self-energies 
and a parton spectral function. 

The PHSD approach incorporates the latter medium 
effects in line with a 1QCD equation-of-state as discussed 
in Section II and also includes a dynamical hadronization 
scheme based on covariant transition rates. As has been 
shown in our previous study [31], the elliptic flow V2 from 
PHSD (red solid lines in Figs. 2 and 3) agrees with the 
data from the STAR and PHENIX Collaborations and 
clearly shows an increase with bombarding energy. As 
was demonstrated in the thorough analysis of Rcf. [27], 
the difference between STAR/PHENIX vi results is less 
than 2%-5% below pt — 2.5 GeV/c. At higher transverse 
momentum, the STAR elliptic flow V2 is systematically 
larger than the PHENIX V2 and the ratio tends to grow 
with pt, reaching the value of 20% at pt — 5.5 GeV/c. 
The differences in V2 at higher pt might be attributed 
to non-flow effects due to di-jets, which are mostly sup- 
pressed by the rapidity gaps in the case of the PHENIX 
measurements. Anyhow, we do not consider such high 
transverse momenta. 

Note that PHSD and AMPT-SM practically give the 
same elliptic flow at the top RHIC energy of ^/snn = 
200 GcV. However, PHSD is more elaborated and in- 
cludes more realistic properties of dynamical quasiparti- 
cles, especially in the vicinity of the critical energy den- 
sity. 

An explanation for the increase in t>2 with collision en- 
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FIG. 3: (Color online) Average elliptic flow V2 of charged 
particles at mid-pseudorapidity for two centrality selections 
calculated within the PHSD (solid curves) and HSD (dashed 
lines) models. The V2 STAR data (stars) for minimal bias 
are the same as in Fig. 2 (stars); the preliminary PHENIX 
data [46] are plotted by filled circles and other data are taken 
from the compilation in Ref. [47]. 



crgy is provided in Fig. 4. Here we show the partonic frac- 
tion of the energy density with respect to the total energy 
where the energy densities are calculated at midrapidity. 
As discussed above, the main contribution to the elliptic 
flow is coming from an initial partonic stage at high y/s. 
The fusion of partons to hadrons or, inversely, the melt- 
ing of hadrons to partonic quasiparticles occurs when the 
local energy density is about e ps 0.5 GeV/fm 3 . As fol- 
lows from Fig. 4, the parton fraction of the total energy 
goes down substantially with decreasing bombarding en- 
ergy while the duration of the partonic phase is roughly 
the same. The maximal fraction reached is the same in 
central and peripheral collisions but the parton evolution 
time is shorter in peripheral collisions. One should recall 
again the important role of the repulsive mean-field po- 
tential for partons in the PHSD model (see Fig. 1) that 
leads to an increase of the flow V2 with respect to HSD 
predictions (cf. also Ref. [48]). We point out in addi- 
tion that the increase of V2 in PHSD relative to HSD is 
also partly due to the higher interaction rates in the par- 
tonic medium because of a lower ratio of r\j s for partonic 
degrees of freedom at energy densities above the criti- 
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FIG. 4: (Color online) Evolution of the parton fraction of 
the total energy density at midrapidity for different collision 
energies at impact parameters 6=1 and 10 fm. 



cal energy density than for hadronic media below the 
critical energy density [49, 50]. The relative increase in 
t>3 and «4 in PHSD essentially is due to the higher par- 
tonic interaction rate and thus to a lower ratio rj/s in the 
partonic medium, which is mandatory to convert initial 
spacial anisotropics to final anisotropies in momentum 
space [51]. 



B. Higher-order flow harmonics 

Depending on the location of the participant nucleons 
in the nucleus at the time of the collision, the actual 
shape of the overlap area may vary: the orientation and 
eccentricity of the ellipse defined by the participants fluc- 
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FIG. 5: (Color online) Projection of a single peripheral 
Au+Au (200 GeV) collision on the transverse plane. Specta- 
tor and participant nucleons are plotted by empty and filled 
circles, respectively. The participant plane transverse axes 
are marked by stars (x*,y*). 
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FIG. 6: (Color online) Average anisotropic flows V3 and v<t 
of charged particles at mid-pseudorapidity for minimum-bias 
Au + Au collisions calculated within the PHSD (solid line) 
and HSD (dashed line) models. 



tuates from event to event. As seen from Fig. 5, due 
to fluctuations the overlap area in a single event can 
have, for example, a rotated triangular rather than an 
almond shape. Note, however, that by averaging over 
many events an almond shape is regained for the same 
impact parameter. 

Recent studies suggest that fluctuations in the initial 
state geometry can generate higher-order flow compo- 
nents [10, 19, 20, 22]. The azimuthal momentum distri- 
bution of the emitted particles is commonly expressed in 
the form of Fourier series as 



E- 



d 3 N 



d 3 p 
d 2 N 



2irpTdpTdy 



1 + V 2v n (p T ) cos[n(V> - *„)] , (4) 



n=X 



where v n is the magnitude of the nth order harmonic 
term relative to the angle of the initial-state spatial plane 
of symmetry The anisotropy in the azimuthal an- 
gle tp is usually characterized by the even-order Fourier 
coefficients with the reaction plane \&„ = f jjp: v n = 
(exp(in(ip — ^>fip))) ( n = 2,4, ...), since for a smooth an- 
gular profile the odd harmonics vanish. For the odd com- 
ponents, say W3, one should take into account event- by- 
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FIG. 7: (Color online) Impact parameter dependence of 
anisotropic flows of charged particles at mid-pseudorapidity 
for minimum-bias collisions of Au+Au at ^snn = 200 GeV. 
Experimental points are from Ref. [55]. 



event fluctuations with respect to the participant plane 
= $pp. We calculate the v 3 coefficients with re- 
spect to * 3 as v 3 {* 3 } = (cos(3[V> - * 3 ]))/Res(* 3 ). The 
event plane angle ^3 and its resolution Res(^I' 3 ) are cal- 
culated as described in Rcf. [52] via the two-sub-events 
method [53, 54]. 

In Fig. 6 we display the PHSD and HSD results for 
the anisotropic flows w 3 and V4 of charged particles at 
mid-pseudorapidity for Au+Au collisions as a function 
of \/snn- The pure hadronic model HSD gives w 3 ~ 
for all energies. Accordingly, the results from PHSD 
(dashed red line) are systematically larger than from 
HSD (dashed blue line). Unfortunately, our statistics are 
not good enough to allow for more precise conclusions. 
The hexadecupole flow V4 stays almost constant in the 
energy range ^snn ~ 10 GeV; at the same time the 
PHSD gives noticeably higher values than HSD, which 
we attribute to the higher interaction rate in the par- 
tonic phase, i.e., a lower ratio of 77/s for the partonic 
degrees of freedom [49, 50]. 

Alongside with the integrated flow coefficients v n the 
PHSD model reasonably describes their distribution over 
centrality or impact parameter b. A specific comparison 
at y/SNjy = 200 GeV is shown in Fig. 7 for V2, w 3 , and 
W4. While «2 increases strongly with b up to peripheral 
collisions, w 3 and v± are only weakly sensitive to the im- 
pact parameter. The triangular flow is always somewhat 
higher than the hexadecupole flow in the whole range of 
impact parameters b. 

Recently, the triangular flow at ^snn = 200 GeV has 
been recalculated in the updated AMPT model [56] . The 
values of model parameters for the Lund string fragmen- 
tation and the parton scattering cross section from the 
previous default version have been refitted to describe 
the charged multiplicity distribution, transverse momen- 
tum spectra, and elliptic flow for Au+Au collisions at 



FIG. 8: (Color online) Time evolution of v n for Au+Au colli- 
sions at y/SNN = 200 GeV with impact parameter 6=8 fm. 



yfsWN = 200 GeV. In the novel AMPT version the par- 
ton scattering cross sections decrease from about 10 to 
1.5 mb. As compared to the old AMPT result w 3 «0.4, 
the new value V3 «0.2 is consistent with the PHSD results 
in Fig. 6. Note, that the magnitude of the PHSD triangu- 
lar flow at y/SNN = 200 GeV is similar also to that from 
the (3 + 1)D viscous hydrodynamical model [57] with the 
specific viscosity r\j s =0.08. In our calculations the low 
transverse momentum particles with px < 1 GeV/c are 
dominating. Unfortunately, experimental data for this 
momentum range are not available. 

Figure 8 shows the time evolution of flow coefficients 
«2, i> 3 , and U4 for a Au+Au collision at impact parameter 
b=8 fm. They reach their asymptotic values by the time 
of 6-8 fm/c after the beginning of the collision, which cor- 
responds to the dominantly partonic phase (cf. Fig. 4). 
Thus, collective flows are formed in the early partonic 
stage of the collision. 

Different harmonics can be related to each other. In 
particular, hydrodynamics predicts that V4 oc (V2) 2 [58]. 
The simplest prediction that 1)4 = 0.5(«2) 2 is given for a 
boosted thermal freeze-out distribution of an ideal fluid, 
in Ref. [16]. In this work it was noted also that 1)4 is 
largely generated by an intrinsic elliptic flow (at least at 
highly) rather than the fourth order moment of the fluid 
flow. This is a motivation for studying the ratio i»4/(f2) 2 
rather than W4 alone. As is seen in Fig. 9, indeed the 
ratio calculated within the PHSD model is practically 
constant in the whole range of ^snn considered, but 
significantly deviates from the ideal-fluid estimate of 0.5. 
This result is qualitatively consistent with the behavior of 
these harmonics in Figs. 3 and 6. In contrast, neglecting 
dynamical quark-gluon degrees of freedom in the HSD 
model, we obtain a monotonous growth of this ratio. 

The dependence of the «4/(i>2) 2 ratio versus the num- 
ber of participants N part is shown in Fig. 10 for charged 
particles produced in Au + Au collisions at ^snn = 
200 GeV. The PHSD results arc roughly in agreement 
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FIG. 9: (Color online) Beam-energy dependence of the ratio 
1)4/(^2) 2 for Au+Au collisions. The solid and dashed curves 
are calculated within the PHSD and HSD models, respec- 
tively. 



FIG. 11: (Color online) Scaling of ^2/^2 as a function of the 
number of participants for different beam energies. Exper- 
imental points are for Au + Au collisions at ^/snn = 200 
(circles) and 62 GeV (squares) [60]. 
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FIG. 10: (Color online) Participant number dependence of 
the V4/(v2) 2 ratio of charged particles for Au+Au ( v / sjvjv = 
200 GeV) collisions. The experimental data points for 0.5< 
Pt <0.8 GeV/c are from Ref. [46]. 



with the experimental data points from Ref. [59] but over- 
shoot them for N part > 250. We will come back to this 
quantity in the last subsection when discussing its px 
dependence. 



C. Scaling of flow coefficients 

The i>2 coefficient measures the response of the heated 
and compressed matter to the spatial deformation in the 
overlap region of colliding nuclei, which is usually quan- 



tified by the eccentricity €2 =< y~ 



> I < x- + y 2 >. 



Since the flow response (1)2) is proportional to the driv- 
ing force (£2), the ratio ^2/62 is used to compare different 
impact parameters and nuclei. 



In Fig. 11 this ratio is plotted as a function of the 
participant multiplicity N part . Note that in these calcu- 
lations the same eccentricity £2 was used as in the exper- 
iment [61]. All PHSD results for ^/s^v ~ 40 GeV are 
very close to each other and in agreement with experi- 
ment. At lower collision energies this scaling starts to be 
violated with decreasing ^/s^n. 

A remarkable property - universal scaling - has been 
proposed in Ref. [62] (see Fig. 12). It appears that 1)2/ &z 
plotted versus (1/ ' S)dN c h/dy falls on a "universal" curve, 
which links very different regimes, ranging from Alternat- 
ing Gradient Synchrotron (AGS) to RHIC energies. Here 
S = it y/ < x 2 >< y 2 > is the overlap area of the collision 
system and dN c h/dy is the rapidity density of charged 
particles. 

As seen in Fig. 12 very different systems vary essen- 
tially only in a single scale; this scale does not depend on 
the collision energy (note that the NA49 data - filled tri- 
angles - are beyond this systematics) and is connected to 
the total entropy produced [64]. Indeed, peripheral and 
low-energy collisions are likely to produce systems with 
incomplete thermalization. Since rescattering of the par- 
ticles is rare in the low-density regime, little change oc- 
curs, on average, to the initial momentum distributions. 
The measured elliptic flow V2 is therefore proportional to 
the initial state eccentricity £2- This quantity and the 
space density of the initial particles dN c h/dy are deter- 
mined via Glauber calculations. In our case we calculate 
these quantities directly within the PHSD model. Thus, 



1 dN ch 

"2 OC —— £ 2 

S dy 



(5) 



We point out that only when event-by-event fluctuations 
in eccentricity are taken into account the universal scal- 
ing is observed in PHSD. The term (1/ S)(dN c f l /dy) con- 
tains information about both the longitudinal structure 
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FIG. 12: (Color online) Scaling of v 2 /e 2 vs (1/ ' S){dN ch /dy). 
The PHSD results are given by lines with open symbols. Pre- 
dictions of ideal boost-invariant hydrodynamics are shown in 
the top panel (from Ref. [63]) and explained in the text. Our 
PHSD results are presented in the bottom panel. The experi- 
mental data points for Au+Au collisions at 200 GeV (circles) 
and 62 GeV (squares) are from Refs. [60, 63]. 



at freeze-out and the final particle number density (which 
is a function of the initial temperature T and baryon 
chemical potential /jb in the hydro limit). 

We therefore use dN c h/dy or, cquivalently, the initial 
entropy density sq in central Au+Au collisions as a proxy 
for the collision energy: At any given collision energy, a 
measurement of dN^/dy in the most central collision 
events fixes the value of so to be used in ideal fluid simu- 
lations at that energy. Assuming linear longitudinal ex- 
pansion without transverse flow at very early time To , the 
quantities dN c h/dy and so are thus related by 



dNch 
dy 



OC T s 



(G) 



In an ideal (isentropic) expansion, the final entropy is 
equal to the initial entropy content of the system [~ the 
initial particle density n(T, /j,b)]- Thus, the systems from 
AGS to RHIC appear to be controlled by a common scale, 
related to the total multiplicity, which varies smoothly 
and drives both v 2 /e2 an d {1 / S){dN / dy) . This conclu- 
sion is a strong indication that microscopic properties of 
the system (equation of state and mean free path) are 



basically unchanged, up to a shift related to this scale, 
in the experimentally addressed energy range. 

In the hydrodynamic limit - implying complete ther- 
malization of the system - the ratio of elliptic flow to ec- 
centricity is saturated at very low impact parameter. In 
this regime the centrality dependence of the elliptic flow 
is mainly determined by the initial elliptic anisotropy of 
the overlap zone in the transverse plane, and the ratio of 
these two should be approximately constant. This is seen 
in Fig. 12 where this correlation is plotted by three hor- 
izontal lines for three different beam energies according 
to Ref. [65]. 

Predictions of ideal boost-invariant hydrodynamics 
based on calculations of Ref. [65] are also presented in 
Fig. 12 (top panel). The lines shown are hydrodynamic 
results for two boost-invariant lattice-inspired equations 
of state (with a quark-hadron phase transition, marked 
as "Q" , and for a pure hadronic system "H" ) calculated 
for the fixed impact parameter (6 = 7 fm) and differ- 
ent particle densities. Note that hydrodynamic results 
do not scale perfectly in this case and in general exhibit 
a somewhat flatter centrality dependence at each colli- 
sion energy. The deviation from experiment is very large 
for peripheral collisions [(1/ ' S)dN c h/ dy ^ 15], where the 
application of hydrodynamics is questionable since the 
mean free path of the degrees of freedom is no longer 
small compared to the transverse size of the system. 

The universal scaling has been investigated in more 
elaborated dissipative hydrodynamic models in Refs. [66, 
67]. A finite shear viscosity rj strongly suppresses the 
buildup of momentum anisotropy and elliptic flow, espe- 
cially for low multiplicity densities. The viscosity effect 
changes the slope of the multiplicity scaling for 1)2/ '£2 but 
preserves, to a good approximation, its general scaling 
with (l/S)dN c h/dy. However, contrary to experiments, 
in hydrodynamic simulations - irrespective of fluid vis- 
cosity - the elliptic flow does not follow the universal 
scaling. In principle, the shear viscosity scaled with the 
entropy density 77/s can be extracted in such an analysis; 
however, the present accuracy of the extracted values is 
not high enough. 

Thus, the experimentally observed scaling in Fig. 12 
puts very strong constraints on the initial microscopic 
properties (entropy density, mean free path, etc.), as well 
as the global longitudinal structure [68]. 



D. Scaling in pseudorapidity 

Another interesting insight on scale invariance in ex- 
perimental observablcs is the longitudinal scaling seen 
experimentally. The multiplicity longitudinal scaling, 
dN c h/dri, is found for a variety of colliding systems and 
is often denoted as "limiting fragmentation". Proposed 
more than 40 years ago by Feynman [69] and Hage- 
dorn [70], this hypothesis implies that the multiplicity 
distribution of particles becomes independent of y/s for 
y/s — > 00. From the microscopic point of view the mul- 



10 




FIG. 13: (Color online) The dependence of the charged par- 
ticle multiplicity on the shifted pseudorapidity. Experimental 
data points are from Ref. [71]. 



tiplicity longitudinal scaling can be understood if the ra- 
pidity distributions of produced hadrons are functions 
of the fraction of the hadron longitudinal momentum 
x = 2pr/y/s alone but not of the total energy. This 
picture is very close to the Bjorken scaling of parton dis- 
tributions. It was found that models combining ideal 
hydrodynamics and hadronic cascades reproduce the lon- 
gitudinal multiplicity scaling pretty well, being rather in- 
sensitive to the "phase" of the system at thermalization. 
This is illustrated in Fig. 13 for very central Au+Au col- 
lisions within the PHSD model. The limiting fragmen- 
tation region is nicely reproduced for (77 — ybeam) ~ — 2, 
while some deviations are seen closer to midrapidity for 
higher collision energies. The situation is different for the 
elliptic flow, since, unlike dNch/drj, the collective flow vi 
is sensitive to the phase of the system as shown before. 

As follows from Fig. 14 the PHSD model reproduces 
the longitudinal V2 scaling up to comparatively low col- 
lision energies (not yet measured). Hadronic results ob- 
tained within in the UrQMD model presented in Ref. [73] 
are quite close to our findings, to be formulated as "a 
qualitative agreement with experiment" . Such a scaling 
was observed also for partonic models such as the AMPT 
and AMPT-SM [73]. 

We mention that the V2 scaling with shifted pseudo- 
rapidity can be obtained also by solving the Boltzmann 
equation with an ellipsoidal profile in the initial trans- 
verse density [74], 



V2 < <Jv > dN c h 
e S dy 



(7) 



where a is the interaction cross section. In line with 
Eq. (5) in the discussion of the universal scaling, such 
an ansatz will naturally lead to an observed-like scaling 
provided that < av > docs not vary with rapidity. 

Furthermore, a comprehensive analysis of the longi- 
tudinal scaling has been performed in Ref. [75] within 
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FIG. 14: (Color online) The V2 dependence versus pseudora- 
pidity as calculated in the PHSD model. The compilation of 
experimental data points is from Ref. [72]. 



simple phenomcnological models, trying various assump- 
tions for hydrodynamic and kinetic descriptions. The au- 
thors conclude that the experimentally observed scaling 
of multiplicity with rapidity and collision energy follows 
from reasonable models of partonic dynamics. Neither 
the free-streaming limit nor the ideal-fluid limit are ex- 
pected to break up this multiplicity scaling [75]. The 
situation is, however, different with the scaling observed 
for the elliptic flow w 2 . It is not clear how this scaling 
could arise within nonideal hydrodynamics, even if its 
initial conditions mirror closely the ones that reproduce 
the scaling observed in dN^/drj. These remarks address 
the shape of the scaling distribution. A more serious 
problem is the absolute value of V2- In terms of Eq. (7), 
to get a reasonable magnitude of V2, the cross section 
a should be increased to the point where the Knudsen 
number is well below unity [75]. As demonstrated above, 
the PHSD model allows us to get reasonable results for 
the multiplicity and V2 longitudinal scaling by default 
without any tuning of parameters. 



E. Differential distributions 

It was shown in the RHIC experiments that, for a 
given centrality, the differential elliptic flow for all ob- 
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FIG. 15: (Color online) Normalized pseudo-rapidity distribu- FIG. 16: (Color online) Directed flow distributions at different 

tions of the directed flow v\ in the transient collision energy centrality for Au + Au collisions at ^Jsnn = 39 GeV. The 

range. The experimental data points are from the STAR Col- experimental data points are from Refs. [78, 83]. 
laboration [78]. 



served hadrons scales to a single curve when plotted as 
vijn q versus Exx/nq, where n q is the number of con- 
stituent quarks in a given hadron species and Etk is 
the transverse kinetic energy for these hadrons [76, 77]. 
This quark number scaling is consistent with the recom- 
bination model, which assumes the collective flow to de- 
velop on the quark level in the QGP phase. Such a scal- 
ing of the elliptic flow z; 2 for identified hadrons has been 
measured by the STAR and PHENIX Collaborations at 
different centralities for top RHIC energies and recently 
also within the BES program [28] . Our paper deals only 
with charged-partiele observables and therefore the num- 
ber of constituent quarks n q is not accurately defined. 
The analysis of scaling properties of identified hadrons 
we postpone to a future study. 

We start with the consideration of rapidity distri- 
butions dv\jdr\ of the directed flow [78] and its beam 
energy dependence as presented in Fig. 15. The di- 
rected flow v\ is the first harmonic coefficient of the 
above Fourier expansion of the final momentum-space 
azimuthal anisotropy Eq. (4), and it reflects the collec- 
tive sidewards motion or "bounce-off" of the particles in 
the final state. Being generated essentially during the 
nuclear passage time ~ 2i?/7, the directed flow probes 
the very early stage of the collision dynamics. In the 
region closer to the beam/target rapidity than to midra- 
pidity, the directed flow is generated very early even at 
a pre-equilibrium stage of the collision [79] and thus it 
probes the onset of bulk collective behavior. In Fig. 15 
the directed flow of charged particles is plotted versus the 
normalized pseudorapidity rj/yteam in the large range of 
the BES collision energies for centrality 30%-60%. We 
observe that vi(r]/yi, eam ) shows a beam-energy scaling 
behavior, though not perfect. Both hydrodynamic and 
nuclear transport models indicate that the directed flow 
is a sensitive signature for a possible phase transition, 



especially in the central region of beam-energies under 
investigation. In particular, the shape of v\(y) in the 
midrapidity region is of special interest because it has 
been argued that differential directed flow may exhibit 
flatness at midrapidity due to a strong, tilted expansion 
of the source. Such tilted expansion gives rise to anti- 
flow [80] . The antiflow is in the opposite direction to the 
repulsive bounce-off motion of nucleons. If the tilted ex- 
pansion is strong enough, it can cancel and even reverse 
the motion in the bounce-off direction and result in a neg- 
ative v\{y) slope at midrapidity, potentially producing a 
wiggle-like structure in v\{y). A wiggle for baryons is a 
possible signature of a phase transition between hadronic 
matter and quark gluon plasma (QGP), although a QGP 
is not the only possible explanation [80-82]. As seen 
from Fig. 15 the slope of the vi(r//yb eam ) distribution at 
r/ = is negative and stays almost constant for ^Jsnn ~ 
10 GeV; its magnitude slightly increases with decreasing 
beam enery, however, exhibiting no irregularities. 

The slope of the pseudorapidity distributions is slightly 
changed when different criteria for centrality selection are 
applied as demonstrated in Fig. 16. The influence of this 
selection is very moderate at midrapidity but becomes 
noticeably stronger in the target-projectile fragmentation 
region with increasing impact parameter. 

Let us continue with differential distributions of the 
elliptic flow f 2 by comparing the pt dependence from 
data with those from the PHSD model. The results from 
PHSD for t>2 (pt) are displayed in Fig. 17 for tJsnn from 
5 to 200 GeV. Also shown are the corresponding results 
from the STAR Collaboration at ^snn = 9, 62, and 
200 GeV (by symbols). The data from PHENIX and 
STAR at midrapidity indicate that the magnitudes and 
trends of the differential elliptic flow [v 2 (pt), centrality 
dependence], are changed only very little over the col- 
lision energy range ^/snn — 62 - 200 GeV, indicating 
an approximate saturation of the excitation function for 
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FIG. 17: (Color online) Beam-energy evolution of transverse 
momentum distributions of Vi(jpt) for Au+Au collisions in 
comparison to the data of the STAR Collaboration from 
Refs. [47, 78, 84]. 



V2 at these energies [27] as exemplified in Fig. 17. We 
mention that the PHSD results underestimate the data 
systematically for pt > 1 GeV, which is attributed to an 
overestimation of scattering of partons with high trans- 
verse momenta. However, the collective flow vi of the 
"bulk matter" is rather well described at all energies 
without any tuning of parameters. 

The momentum distribution of the flow V3 is also 
in a quite reasonable agreement with experiment (see 
Fig. 18). 

As pointed out before, the ratios of flow coefficients 
might shed valuable light on the actual dynamics, since 
especially the ratio V4/(v2) 2 is sensitive to the micro- 
scopic dynamics. In this respect we show the transverse 
momentum dependence of the ratio V4/(v2) 2 in Fig. 19 
for charged particles produced in Au+Au collisions at 
^/sWn = 200 GeV (20-30% centrality). The PHSD re- 
sults are quite close to the experimental data points from 
Ref. [59]; however, they overestimate the measurements 
by up to 20%. The hydrodynamic results - plotted in 
the same figure - significantly underestimate the exper- 
imental data and noticeably depend on viscosity. The 
partonic AMPT model [86] discussed above also predicts 
a slightly lower ratio than the measured one; however, it 
is in agreement with both hydrodynamic models for px ~ 
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FIG. 18: (Color online) Triangular flow U3 as a function of 
transverse momentum pr for Au + Au collision at ^snn = 
200 GeV. The data points are from Ref. [52]. 



2.5 



2.0 



1.5 



> 1.0 



0.5 



0.0 



Au + Au, 200 GeV, 20-30% 
all charged, lr|l < 1.3 



STAR preliminary 
AMPT L.Chen et al. 

Hydro P.Kolb et al. 

Hydro N.Borghini et al. 

— e— PHSD 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
P T [GeV/c] 



FIG. 19: (Color online) Transverse momentum dependence of 
the ratio V4/(v2) 2 of charged particles for Au+Au (at ^/snn = 
200 GeV) collisions. The dashed and dot-dashed lines are cal- 
culated within the hydrodynamic approaches from Refs. [58] 
and [85], respectively. The shaded region corresponds to the 
results from the AMPT model [86]. The experimental data 
points are from the STAR Collaboration [59]. 



0.8GcV/c. Our interpretation of Fig. 19 is as follows: the 
data are not compatible with ideal hydrodynamics and a 
finite shear viscosity is mandatory (in viscous hydrody- 
namics) to come closer to the experimental observations. 
The kinetic approaches AMPT and PHSD perform better 
but either overestimate (in AMPT) or slightly underesti- 
mate (in PHSD) the scattering rate of soft particles. An 
explicit study of the centrality dependence of these ratios 
should provide further valuable information. 
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IV. CONCLUSIONS 

In summary, relativistic collisions of Au+Au from 
^s NN — 5 to 200 GcV have been studied within the 
PHSD approach, which includes the dynamics of explicit 
partonic degrees of freedom as well as dynamical local 
transition rates from partons to hadrons and also the fi- 
nal hadronic scatterings. Whereas earlier studies have 
been carried out for longitudinal rapidity distributions 
of various hadrons, their transverse mass spectra, and 
the elliptic flow V2 as compared to available data at SPS 
and RHIC energies [32, 33], here we have focused on the 
PHSD results for the collective flow coefficients v±, V2, 
i>3, and t>4 in comparison to recent experimental data 
in the large energy range from the RHIC Beam Energy 
Scan (BES) program, as well as different theoretical ap- 
proaches ranging from hadronic transport models to ideal 
and viscous hydrodynamics. We mention explicitly that 
the PHSD model from Ref. [33] has been used for all 
calculations performed in this study and no tuning (or 
change) of model parameters has been performed. 

We have found that the anisotropic flows - elliptic V2, 
triangular and hcxadecapole i>4 - are reasonably de- 
scribed within the PHSD model in the whole transient 
energy range, naturally connecting the hadronic pro- 
cesses at lower energies with ultrarelativistic collisions 
where the quark-gluon degrees of freedom become domi- 
nant. The smooth growth of the elliptic flow t>2 with the 
collision energy demonstrates the increasing importance 
of partonic degrees of freedom. This feature is repro- 
duced by neither hadron-string based kinetic models nor 
a multi phase transport (AMPT) model treating the par- 
tonic phase in a simplified manner. Other signatures of 
the transverse collective flow, the higher-order harmon- 
ics of the transverse anisotropy V3 and U4, change only 
weakly from ^s^n ~ 7 GcV to the top RHIC energy 
of yfsNN = 200 GcV, roughly in agreement with experi- 
ment. As shown in this study, this success is related to a 
consistent treatment of the interacting partonic phase in 
PHSD, whose fraction increases with the collision energy. 

The observables calculated within the PHSD model ex- 
hibit some scaling properties for collision energies above 
\J snn = 40 GeV. In particular, the universal scaling 
of V2/t2 versus (1/ ' S)dN c h/dy (cf. Fig. 12) is approxi- 
mately reproduced as well as the longitudinal scaling of 
the charged particle pseudorapidity distributions of the 
elliptic flow i>2 in (rj — ybeam) representation (cf. Figs. 13 
and 14) in this energy range. This feature is not repro- 
duced by hadronic transport models (such as HSD and 
UrQMD) and meets (severe) problems in the various hy- 
drodynamic descriptions. 

The analysis of correlations between particles emitted 
in ultrarelativistic heavy-ion collisions at large relative 
rapidity has revealed an azimuthal structure that can be 
interpreted as being solely due to collective flow [87. 88]. 
This interesting new phenomenon, denoted as triangu- 



lar flow, results from initial-state fluctuations and a sub- 
sequent hydrodynamic-like evolution. Unlike the usual 
directed flow, this phenomenon has no correlation with 
the reaction plane and should depend weakly on rapid- 
ity. Event-by-event hydrodynamics [89] has been a natu- 
ral framework for studying this triangular collective flow 
but it has been of interest also to investigate these corre- 
lations in terms of the PHSD model. We have found the 
third harmonics to increase steadily in PHSD with bom- 
barding energy. The coefficient ^3 is compatible with zero 
for y/sNN > 20 GeV in case of the hadronic transport 
model HSD, which does not develop "ridge-like" correla- 
tions. In this energy range PHSD gives a positive i>3 due 
to dominant partonic interactions. 

Different harmonics can be related to each other 
and, in particular, hydrodynamics predicts that V4 oc 
(1*2) 2 [58]. In this work it was noted also that V4 is largely 
generated by an intrinsic elliptic flow (at least at high px) 
rather than the fourth-order moment of the fluid flow. 
Indeed, the ratio v^/fa) 2 calculated within the PHSD 
model is approximately constant in the whole considered 
range of ^/snn but significantly deviates from the ideal 
fluid estimate of 0.5. In contrast, neglecting dynamical 
quark-gluon degrees of freedom in the HSD model, we 
obtain a monotonous growth of this ratio. 

The transverse momentum dependence of the ratio 
vaI(v2) 2 at the top RHIC energy has given further inter- 
esting information (cf. Fig. 19) by comparing the various 
model results to the data from STAR, which are inter- 
preted as follows: the STAR data are not compatible 
with ideal hydrodynamics and a finite shear viscosity is 
mandatory (in viscous hydrodynamics) to come closer ot 
the experimental ratio observed. The kinetic approaches 
AMPT and PHSD perform better but either overestimate 
(in AMPT) or slightly underestimate the scattering rate 
of soft particles (in PHSD). An explicit study of the cen- 
trality dependence of these ratios should provide further 
valuable information. 

It will be promising to extend our studies to asym- 
metric heavy-ion collisions that can be used to constrain 
models dealing with flow fluctuations in heavy-ion colli- 
sions but with a larger sensitivity for ^-related observ- 
ables than for V3 [90]. Independently, an extension of 
the PHSD approach to LHC energies with possible color- 
glass-condcnsate initial conditions has to be performed 
in future. 
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